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Abstract. Influence of the chiral symmetry on the many-body problem in multilayer graphene 
in magnetic fields is investigated. For a spinless electron model on the honeycomb lattice the 
many-body ground state is shown to be a doubly-degenerate chiral condensate irrespective of the 
number of layers. The energy spectrum calculated numerically with the exact diagonalization 
method reveals for ABC-stacked multilayer graphenes that the many-body gap decreases 
monotonically with the number of layers. 



1. Introduction 

Graphene's intriguing physics [U [21 [3l |3] is intimately related to the chiral symmetry, a 
fundamental feature of the honeycomb lattice. The symmetry is defined by the anticommutation 
relation, {T-L^V} = 0, with the chiral operator V acts as 

VciV = +Ci{-Ci) for ie.(o), r2 = l, (1) 

where • and o respectively denote the two sublattices. The symmetry guarantees [3 E] the 
topological stability of the doubled Dirac cones, an example of Nielsen-Ninomiya theorem [3 [5], 
as well as the degeneracy of the n = Landau level (LL), which accommodates Aharonov- 
Casher's argument [U]. A remarkable phenomenon in the latter is the delta- function-like density 
of states (DOS) for the n = LL even in the presence of ripples with wavelengths exceeding 
a few lattice constants [10] . These properties are expected to be inherited by multilayer 
graphene [HI [121 [E]) since the leading hopping matrix elements, 70,71 and 73, conserves the 




1 2 3 4 5 6 7 
(bottom) layer index 



9 10 
(top) 



Figure 1. Layer-by-layer charge distribution in the chiral 
condensates for an ABC stack of 10 graphene sheets. 
71/70 = 0.5, where only the electrons in the zero-energy 
LL is shown. The filled (empty) circles, triangles and 
squares correspond to for (p = 3/15^,2/15^ and 

1/15^, respectively. 



bipartite lattice structure, and hence the chiral symmetry. Indeed, a recent theory [2] predicts 
that AB-staked bilayer graphene retains the anomaly in DOS, which even extends to the case 
of a perpendicular electric field that in fact breaks the chiral symmetry. 

Then an interesting question is how the chiral symmetry affects the many-body problem in 
multilayer graphene. Many-body effects have been extensively studied [U [TJl HSl [13 ESI El EDI 
ISlEaiMllMlEaiMlEZllMllMllMlEIllMllMllMllMllMlE^ from the beginning of the 

graphene research due to interests that include the gap opening in the central LL [HJ 351 HSl 33] 
and the unconventional insulating 1/ = state j45j . So far, the role of the chiral symmetry in the 
many-body problem has been studied in the context of the lattice gauge theory, where monolayer 
graphene without magnetic fields is discussed, mainly based on the continuum model [4611471 HH]. 
In the present work, we consider multilayer graphenes in magnetic fields, and investigate the 
many-body problem in terms of the chiral symmetry. 

2. Ground states 

We consider spinless electrons in multilayer graphenes in a magnetic field, for which the kinetic 
energy is 

where jij is the electron hopping between sites {ij). For the nearest-neighbor in-plane sites 
lij = 70) while for the inter-layer vertical hopping 'jij = 71. The magnetic field is included as 
the Peierls phase 6ij such that magnetic flux per elementary hexagon equals to 4> = 6ij 
in units of the magnetic flux quantum h/e. In the string gauge [49j . the magnetic flux 
reads (j) = n/Nc with an integer n and the number of unit cells, Nc- Since the 70-71 model 
preserves the chiral symmetry, one can construct a zero-energy multiplet, V = (V'+iV'-) with 
V'lt = (V'i±) ■ ■ ■ ) '4'M±±), where V'mi's denote M-t-fold degenerate eigenstates of F corresponding 
to chirality ±1, i.e., F-^mi = iV'mit- For a A'^^-layer bulk system with a magnetic flux (j) = n/Nc 
(n = 1,2,---), the degeneracy reads = M_ = nNi. It should be noted that the zero 
modes with chirality -|- or — are localized on sublattice • or o, respectively. Especially for ABC- 
stacked multilayers, the zero modes show sublattiece-selected charge accumulation towards the 
top or bottom layer, which decays rapidly away from the surface [50 \ 1511 [52] . The surface states 
are analogous to the zigzag edge states in the ribbon structure [53] and should be observed 
experimentally as local density of states near zero energy with a scanning tunneling microscope. 

For a half- filled multilayer graphene, the interaction for spinless electrons can be written in 
a particle-hole symmetric form as 

n-mt = J2 f'^* ~ ^) f'^^' ~ ^) ^iji4<^]<^jCi + CiCjcyi) + const (3) 
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Figure 2. Energy spectra of AB-stacked bilayer graphene (a), ABC-stacked trilayer graphene 
(b) and ABA-stacked trilayer graphene (c), plotted against the total chirality, x- Results for 
iVc = 152, 4> = 3/152, 71/70 = 0.5, Vi/Vo = 0.5 are shown. 

with interaction strength Vij. Since it is hard to treat all the many-body states exactly, we 
focus on the many-body problem within the zero-energy LL. Such a treatment is adequate as 
long as the interaction energy is sufficiently smaller than the Landau gap between the central 
and the adjacent LLs. The projection onto the zero-energy LL is performed as c = (V'V'^)c, 
where c = (ci, • • • ,cn)'^ with the number of sites A^, and denotes the projection matrix. 
Note that the projected creation and annihilation operators are no longer fermionic; they obey 
anticommutation relations {ci,cj} = (V'V'^)ii and {ci,Cj} = {c|,c^} = 0. In what follows, we 
assume that the electron-electron repulsion only acts between electrons on different sublattices 
[Vij > for = (•, o)]. The simplest example is the nearest-neighbor repulsion, which is the 
leading term for spinless electrons. The total Hamiltonian in the projected subspace reads, up 
to a constant, 

iG'j'Go klmn 

where dl^j^ = c^'ipj^^^ is the creation operator for the zero mode ipm±, while the 
pseudopotential appearing for d,d^ is defined as Vkimn = I]je.,ieo ^u(V'fc)I(V'i)j (V'm)j('0n)^• 
Clearly, H is semi-positive definite; for an arbitrary many-body state |$), (<^|?^|<I>) = 
I I]je.,ieo ^ij IfcjCj + (cjCj)'l"]|^>)| > 0. Since V. commutes with Q = YjnMm+dm.+ - d^rn-dm-), 
eigenstates of 7i are classified according to the total chirality, x = (G)- The highest and lowest 
sectors in x correspond to chiral condensates, 

|G±) = 4± • • • 4±±l^<). S\G±) = ±M±|G±), (5) 

where all the zero modes with chirality + (or — ) are occupied, while others are empty. The 
chiral condensates are in fact the many-body ground states, since they are annihilated by both 
of c?|,_|_(i|-|- and dm^dn±-, hence {G±\ii\G±) = 0. Thus the chiral condensates form a ground-state 
doublet, ^ = (|G+), so that any unitary-mixture, ^ ^ ^^oj with u € U(2), is again an 

eigenstate. Since the above discussion holds irrespective of the number of layers, it is expected 
that ^ has properties similar to those in the monolayer case, such as vanishing Hall conductivity, 
or the Kekule-type bond-order along armchair edges [MIES]- 

As a physical property peculiar to the multilayer case, we plot in Figure [T] layer-by-layer 
charge distribution in \G±) for an ABC-stacked multilayer with N^, = 15^ and Ni = 10, 
where only the electrons in the zero-energy LL is shown. Since an experimental cyclotron 
frequency Uc is too small to treat numerically based on our lattice model, we adopt a rather 




Figure 3. Energy gap A of ABC-stacked multilayer 
graphene composed of Ni graphene sheets as a function of 
Ni. Here a result for 71/70 = 0.5, Vi/Vq = 0.5 is shown. 
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large parameter 71/70 = 0.5 to access the realistic parameter region 71 > Uc- One can see charge 
accumulation towards the top or bottom layer, which reflects the property of the zero modes 
as the surface states [SUl ED E2]- Notably, a shoulder structure evolves on the intermediate 
layers with increasing (f>, Although the non-monotonic behavior disagrees with the results of the 
continuum model, where low-energy wave functions decays exponentially away from a surface, 
the discrepancy seems to be due to the finite size effect of our lattice model. 

3. Excited states 

Next we investigate excited states at half filling by numerically diagonalizing the projected 
Hamiltonian (jH). For simplicity, only the nearest-neighbor repulsions in the same layer Vq and 
those between adjacent layers Vi are taken into account. Fig. [2]shows typical examples of energy 
spectra plotted against the total chirality, Xi for AB-stacked bilayer graphene and for ABC- or 
ABA-stacked trilayer graphene. The spectra are symmetric with respect to X) since now we have 
the degeneracies of the zero modes with chirality -|- and — are equal. From this total-chirality- 
resolved plot one can confirm that the ground states are indeed chiral condensates, as expected 
from the above discussion. More importantly, these data suggest that the first excited states 
are obtained by flipping a single-chirality of the chiral condensates as 

\E±)=Y,Cmndl^dn±\G±) (6) 

mn 

with coefficients Cmn determined numerically. Although the size of the Hilbert space blows up 
exponentially with the number of layers, the ansatz ([6]) enables us to calculate the energy gap A 
between the ground state and the first-excited state by restricting ourselves to the sector of total 
chirality x = i(-^± ~ 2). Figure [3] plots thus obtained energy gap A of ABC-stacked multilayer 
graphene composed of Ni graphene sheets. The decrease in A as a function of Ni is consistent 
with the fact the overlap of the opposite-chirality zero modes, which are accumulated near the 
opposite surfaces, becomes small with increasing Ni. The saturation in A for large Ni suggests 
that the overlap does not completely vanish due to the finite size effect. 

4. Summary 

We have considered spinless fermions in multilayer graphene and investigated many-body states 
in the quantum Hall regime. For the electron-electron repulsion restricted to unlike sublattices, 
the ground state is exactly identified to be chiral condensates as in the case of monolayer 
graphene. Using the exact diagonalization method, we have numerically calculated the many- 
body gap, which decreases monotonically with increasing the number of layers. 
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